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Abstract 

■ A first-principles derivation is given of the imaginary part of a Green's function in real-time thermal 
field theory. The analysis and its conclusions are simpler than in the usual circled-vertex formalism. 
The relationship to Cutkosky-like cutting rules is explained. 

< 

In a recent paper, Bedaque, Das and NaikW have discussed some cancellations among diagrams that 
{Sj ■ contribute to the imaginary parts of finite-temperature amplitudes. A purpose of this note is to point 
out that any such cancellations occur because the usual treatment of the imaginary parts, which 
starts in x space and introduces circled vertices^, is unnecessarily complicated. A momentum-space 



(N 



treatment, which starts from first principles^] , and stays close to the physics leads to the answer much 



O 



43 
Oh 



more simply and directly. 



Suppose, for definiteness, that we want to calculated dilepton production from a plasma at temper- 
ature 1/(3. The plasma is defined to be composed of quarks and gluons only, with any leptons or 



photons that may be produced immediately escaping from it. The initial density matrix is 



Po = Z^VeM-pH) (la) 

where H is the hadronic part of the Heisenberg-picture Hamiltonian and V is a projection operator 
which removes any states that do not contain just quarks and gluons. We may choose to express V in 
^ ' terms of a complete orthonormal set of plasma states \i in), so that 



PQ 



Z- 1 2^|iin)(iin|e" /SJ;r (16) 



The partition function Z is defined in the usual way, so as to make trpo = 1- The probability of 
emission of a lepton pair of total momentum q is then calculated from 

W^(q) = f d 4 xe^ x ^(f out|J M (x)p ^(0)|/ out) (2) 
J f 

where J(x) is the hadronic part of the Heisenberg-picture electromagnetic current and \f out) are a 
complete set of final plasma states. 

The usual zero-temperature perturbation theory applies to the matrix elements that arise in the sum 
(2). Introduce first an interaction picture that coincides with the Heisenberg picture at time to'. 

J(x) = U t0 (t ,x°) J ta (x)U to (x°,to) (3) 
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where ^ 

U to (t 2 ,h) = A(t 2 )A- 1 (ti) =Texp dt J d 3 xH^ T (t,x 

_ e i(t-t )Hl° e i(t-t )H ^ 

Here T denotes ordinary time-ordering and H^ T is the interaction Hamiltonian density written as a 
functional of interaction-picture fields and canonical momenta. We have 

e -f>H =e -K° U^Qo-ifrto) (5) 

We take the limit to — > — oo, so that the interaction-picture fields become the in fields. We use also 

|/out)=C/ in (-oo,oo)|/in) (6) 

Then 

W^(q) = Z- 1 J d A xe iq - xS ^{f in|[/ in (oo,x°) J™(x) U in (x°, -oo)\i in) 

{i m\e-P H ° n U in (-oo - i0,O) J?(0) f/ in (0,oo)|/ in) (7) 
By using completeness for the states / we obtain: 



Z' 1 J d 4 xe tq - x Y^ 



(i in|e-^» n f/ in (-oo-i/3,0) Jl n (0) U in (0,x°)Jl n (x) U in (x°, -oo)\i in) (8a) 
We recognise thist 5 ! as the standard perturbation-theory expression for 

GfM = J d 4 xe^tr(p J,(0)J M (x)) (9a) 

with the Keldysh contour^], which in the complex t plane runs along the real axis from — oo to +oo, 
back along the real axis to — oo, and then down to — oo — i(3. We may obtain an alternative expression 
by using instead completeness of the states i in (6b), and the fact that H™ is the time-translation 
operator for the in fields. We find, instead of (9a), 

e'^GfM = e~ Pq ° J d 4 xe^-tr(p o J At (x)J,(0)) (96) 

For the conserved hermitean electromagnetic current, G 2 } l ,{q) is symmetric in fi, u and is real. Of 
course, the equality of (9a) and (9b) is nothing but the standard relational 

G 12 (t,x) = G 21 {t-ip,x) (9c) 

and one of its consequences is that 

GfM) = Jtfxe^trfaTJpWMO)) (10a) 

satisfies 



Im^( 5 ) = i(e rf + C( 9 ) (106) 
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We shall find that there are simple generalised Feynman rules to calculate Gj^(q), and that they 
may be recast into a formalism resembling the zero-temperature Cutkosky formula^], though with 
important differences^]. 

To derive these, we start with the expression (8a) for Gll(q). It is well-known! 9 ] that one may replace 
the argument (— oo — i(5) in the first U with — oo, so long as one does the same in the calculation 
of Z. That is, in the original definition (la) of po we may replace H with H . Of course, the full 
Hamiltonian is still used to describe how the system evolves away from its initial density matrix p . 

The physical reason that allows this modification of the initial density matrix is as follows. The photon 
emission rate which we wish to calculate is constant in time and depends only on the density matrix 
at the given time. If we change the density matrix in the remote past, by multiplying the interaction 
by the usual adiabatic switching factor e _e '*' used in ordinary scattering theory, this will not affect 
the photon emission rate at t = 0. In practice, we calculate the total emission over a very long time 
interval T, and then divide by T by dropping a factor 6(0). For most of the time interval T the 
switching factor has little effect, and its presence just changes the result by a term of order eT and so 
it is unimportant. 

With this change, 

Gf v {q) = Z- 1 j d A xe iq - x ^2e-^(i in|[/ in (-oo, 0) J^ n (0) C/ in (0, oo)|/ in) 
J i,f 

{f m\U in (oo, x°) Jl n (x) U in (x , -oo)\i in) (11) 

The double sum (11) may be calculated from generalised Feynman diagrams that use the Keldysh- 
contour 2x2 matrix propagator. In the case of a scalar field this matrix propagator is 

^ - ( 2 ^ w m2) ) + ^ - ™ 2 > »<*■> ( ! ! ) < I2 "» 

where Dp is the zero-temperature Feynman propagator and n is the Bose distribution 

There is a similar matrix propagator for the spin-i case. For gauge fields, we use the formalism^] 
in which only the two physical degrees of freedom of the gauge field have a thermal term; the other 
two, and the Faddeev-Popov ghost field, have only the first matrix in (12a). However, for simplicity 
pretend first that the fields of the quarks and gluons that make up the plasma are scalar. 

To derive the diagrammatic expansion of (11), expand the U's as in (4) and apply the zero-temperature 
Wick theorem to each matrix element in the double sum. Each product of operators UJU in the second 
matrix element is written as a sum of products of zero-temperature Feynman propagators {0\T<p<p\0) 
and normal-ordered creation and annihilation operators. For the first matrix element we instead have 
{0\T<p<p\0) and therefore complex conjugates of propagators. In order to evaluate the matrix elements 
of the normal products, we must introduce the finite volume V of the plasma, and a discrete spectrum 
for the fields: _ 

^{x) = ^2 a r exp(— iuj r x° + ik r .x) + h c (13a) 

r V r 

with 

[a r ,aj] = S rs (136) 

The states are then labelled by the occupation numbers n r of the various single-particle modes r. In 
the V — > oo continuum limit 
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[a(k), at (k')] = (2tt) 3 2fc° <5 (3) (k - k') (14) 
but we do not take this limit until we have used the standard relations 

a r \ni,n 2 • • •) = \/n^\n 1 ,n 2 . . . n r - 1 . . .), aj|ni,n 2 • • •) = \Av + l|ni,n 2 . . . n r + 1 . . .) (15) 

In order to give a non-zero contribution, the indices r on the operators a and aft in (11) must be equal 
in pairs. The possible pairings are 

(i) the index on one of the operators a in the normal product in the second matrix element matches 
that on one of the operators aft in the anti-normal product in the first matrix element 

(n) the index on one of the operators aft in the second matrix element matches that on one of the 
operators a in the first matrix element 

(m) the index on one of the operators a in the first matrix element matches that on one of the 
operators aft in the same matrix element 

(iv) the index on one of the operators a in the second matrix element matches that on one of the 
operators aft in the same matrix element 



We need not consider the case where three or more indices are equal to each other; this would reduce 
the number of indices left to be summed and sot 11 ! , according to (14), give fewer powers of V when 
we go to the continuum limit. (For a system in thermal equilibrium these nonleading powers actually 

cancel^, though this is not true more generally.) For each pair of equal indices r, we use (15) to 
replace the corresponding two operators with n r in the cases (i) ,(iv) , and n r + 1 in case (ii). We 
then sum over n r using 



ne 



) 



n(to) 



(16) 



where n(oj) is again the Bose distribution. 





(a) 



(b) 



Figure 1: (a) a matrix element that contributes to (11), with (b) its square summed over the momenta. The 
lines at the top of the diagrams are the spectator particles in the heat bath. 
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Figure 2: (a) The diagram of figure lb with one of the internal lines k shown explicitly; in (b) and (c) this line 
instead is a heat-bath particle. 

We then take the continuum limit and represent the results diagrammatically. Consider first the 
contribution from cases (i) and (ii) only. In (11) we have a matrix element, shown in figure la, which 
has to be squared and integrations applied to the momenta on the incoming and outgoing lines. The 
lines at the top of the diagram are the spectator particles of the heat bath. In the squared matrix 
element the incoming lines, those on the left in figure la, correspond to case (i); with each of them is 
associated 2ir5 + (k 2 — m 2 )n(k°). The outgoing lines on the right correspond to case (ii); with them is 
associated 2n8 + (k 2 — m 2 )(n(k°) + 1). The result of squaring figure la and performing the momentum 
integrations may be depicted as in figure lb, where 2 labels the first matrix element in (11) and 1 
the second. To agree with what has been said, the lines joining them correspond to the {21} element 
of the thermal propagator matrix iD of (12a); the direction of the arrow denotes whether the energy 
flow from 2 to 1 is positive or negative. The matrix element 1 is an ordinary zero-temperature matrix 
element; its internal lines are zero-temperature Feynman propagators iDp- The matrix element 2 is 
its complex conjugate, so its internal lines are rather (iDp)* ■ 

In figure 2a we expose one such internal line k in the matrix element 1. Figures 2b and 2c show 
diagrams corresponding to case (Hi). The internal line k is replaced with two external lines of the 
same momentum. In the matrix element 1, this corresponds to a heat-bath particle k undergoing 
forward scattering and returning to the heat bath with the same momentum. This matrix element 
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(b) 

Figure 3: (a) an "uncuttable graph", with an example (b) of a pair of interfering physical processes to which it 
corresponds 

is multiplied by a corresponding matrix element 2 in which the heat-bath particle k is a spectator. 
Figure 2b may be obtained from figure 2a by replacing iDp(k) with 2ir5 + (k 2 — m 2 )n(k°), and for 
figure 2c we instead need 2ttS~ (k 2 — m 2 )n(k°). If we add the three figures together, iDp(k) becomes 
just the {11} element of the thermal propagator matrix iD. Exactly similar arguments apply to case 
(iv) and the matrix element 2, where instead we arrive at the {22} element of iD. So we may omit 
graphs where there is interference with spectator particles in the heat bath, and instead use thermal 
propagators on internal lines of the matrix elements 1 and 2. Then we need no longer explicitly draw 
the heat-bath spectators in our diagrams. 

Notice that the matrix elements 1 and 2 need not both be connected. This can lead to the "uncuttable" 
diagrams of Kobes and Semenoff ® . An example of such an uncuttable graph is figure 3a. One of the 
physical processes to which it corresponds is shown in figure 3b, where a connected graph on the left 
interferes with a disconnected one on the right. 

The prescription, then, is that the incoming-current vertex is of type 1 and the outgoing one of type 
2. For the internal vertices we sum over either possibility. Of course, some of the resulting graphs 
vanish for kinematic reasons. An example is figure 4: one of its internal vertices is a type-1 vertex 
connected just to three type-2 vertices. Because the {21} element of iT) puts the corresponding line 
on-shell, and - at least in the case of nonzero mass - it is kinematically impossible to have a three-line 
vertex where all three lines are on shell, the graph vanishes. 

The graphs of figure 5 contain "self-energy" insertions. Figures 5a and 5b are associated with forward 
scattering on a heat-bath particle, as has been explained above, while figure 5c represents either the 
absorption of a particle from the heat bath or the emission of an extra one into it. The graph of figure 
5d vanishes for the kinematic reasons explained above. These self-energy diagrams are delicate, in 
that two singular propagators having the same argument are multiplied together, but by switching 
off the interaction in the remote past and future one may showt^l that one should handle them in a 
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Figure 4: a graph that vanishes for kinematic reasons 



way similar to zero-temperature self-energy insertions. The thermal propagator i~D(k) of (12a) may 
be written asN 

tD(k) = M(k)(iD(k))M(k) 



' /3|fc°|/2 e -(5k°/2 

M(k) = ^n(k ) ( ePk o /2 ^| fe o |/2 



iD(fc) 



iT> F {k) 
—iD* F {k) 



(17a) 



where Dp is, as usual, the zero-temperature Feynman propagator. The matrix self energy II(fc) has 
the form 

U(k) = M" 1 (A;)n(A;)M- 1 (A;) 



lL[k) ~ v o -u*(k) 



(17b) 





Figure 5: graphs with self-energy insertions 



In consequence, the Dyson sum of repeated self-energy insertions gives the dressed thermal propagator 

iD'(fc) = M(k) (it)'(k)) M{k) 



it)'(k) 



i(k 2 -m 2 -n(fc) -Me)" 1 







-i{k 2 -m 2 -U*(k) -ie)' 1 

To calculate the sum of graphs in figure 5 we need for the upper line with its insertions 

1 1 



(17c) 



*[D'(fc)] 21 



k 2 — m 2 — n(/c) + ie k 2 — m 2 — II* (k) — ie 

For values of k for which H(k) is real, this is 

i[B'(k)} 21 = 2ir5(k 2 -m 2 - U{k)) (n(k°) + 9(k )) 
= 2irZp5(k 2 - m 2 p (k )) (n{k°) + 6(k )) 



n(k°) (eP k °e(k°) + #(-fc )) (18a) 



(186) 



where m|(fc°) is the value of k 2 for which (k 2 — m 2 — U(k)) vanishes. The zero-temperature values 
of mp and Zp are just the renormalised mass and the contribution from the self energy to the charge 
renormalisation constant. These zero-temperature values arise just from figures 5a and 5b, not 5c. At 
nonzero temperature the additional contributions to mp and Zp may be thought of as providing extra 
shifts to the mass and the charge. Usually these are fc°-dependent; they still arise only from figures 5a 
and 5b because, as can be checked from (17b), II 21 vanishes when II is real. When H(k) is complex, 
all three diagrams of figure 5 contribute. 

Similar considerations apply to self-energy insertions in the 11 and 22 elements of the thermal propa- 
gator D, which occur within the 1 and 2 matrix elements in figure lb. 

The extension of this treatment to include spin-^ particles is obvious. In the case of gauge particles, 
note that the sum over states in (lb) should extend only over physical states i, since the statistical 
mechanics from which it is derived deals only with the possible physical states of the ensemble. The 
summation (2) over states is again a summation over probabilities of achieving any physical state j. 
However, here we may include also unphysical states, so as to make up a complete set. This is because 
a process that begins from a physical state has total probability zero to end up in an unphysical state: 
in the canonical quantisation the Faddeev-Popov ghosts are introduced expressly to achieve thi S [ 15 L 
Hence we again arrive at the expression (9a), though with the trace understood to be restricted to a 
summation over the physical states i. In a similar way, if we restrict the states j to be physical, we may 
extend the set of states i to make up a complete set, because the ghosts also ensure that a process 
that finishes in a physical state has zero total probability that it began in an unphysical one. So 
again we may derive (9b). So all our diagrammatic analysis is still valid, provided that any external 
gauge particles in the graphs are restricted to their physical degrees of freedom, with no external 
ghost lines, and for internal gauge particles only the physical degrees of freedom are thermalised, with 
their unphysical modes and the ghosts having just the zero-temperature Feynman propagators^ . Of 
course, when zero-mass gauge particles are present there are various infrared divergences. These must 
cancel in physical quantities, and the gauge invariance together with unitarity can also cause other 
cancellations ^ '^l . 
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